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Abstract 

Diffusion of particles in velocity space undergoing turbulent field was extensively stud- 
ied in the problem of warm beam relaxation. Under low field intensities the diffusion is 
described by the Fokker-Planck equation with the diffusion coefficient given by quasilin- 
ear theory. This diffusion coefficient is calculated on the free particle propagator and for 
weak fields its renormalization due to orbit diffusion is not necessary. To study effects 
which should be taken into account when the intensity of the turbulent field is increased a 
numerical simulation of particle motion in the external field of Langmuir waves with given 
/c-spectrum and random phases is done. For strong fields meansquare velocity evolution 
shows that ballistic regime in the very beginning is changed for oscillatory one in the 
intermediate stage and later for diffusion. Asymptotically it behaves like fractional power 
of elapsed time with the exponent dependent on the particular field spectrum. Such evolu- 
tion in the whole temporal interval of simulation is recovered from the numerical solution 
of generalized Fokker-Planck equation with time dependent diffusion coefficient obtained 
from the microscopic approach. The analytical approximation for this solution is also 
given. 

1 Introduction 

Diffusion of particles in random external fields could be considered in relation to a general prob- 
lem of transport in plasmas. The test particle approach, which is simpler and more controllable 
than selfconsistent one, helps to analyze some particular aspects of turbulent transport. 

One of not clarified issues reported in papers [1,2] concerns the enhancement of a diffusion 
coefficient in velocity space as compared to its quasilinear value. It was supposed [1] that the 
enhancement was caused by peaks in spectrum being formed in a course of evolution. The 
tendency of formation of nonuniform structures in plasma-beam system was pointed out in 
the works [3]. Here we are interested in what effects are to be taken into account when the 
spectrum of external field become stronger and narrower, i.e. more peaked. 

We have made a direct simulation of test particle motion in prescribed random fields for 
a variety of spectra, and found from obtained data the evolution of average velocity and dis- 
persion. Then the generalization of the Fokker-Planck equation were considered in order the 
solutions would be consistent with simulation in wide range of variation of external field spectra. 

As it was expected for low intensity and broad spectrum (small Kubo number) solutions of 
Fokker-Planck equation with quasilinear diffusion coefficient agree with the results of simula- 
tion. When the form of spectrum is taken peaked (Kubo number becomes larger than the unit) 
the velocity dispersion grows very fast at a time less than the field correlation time, and such 
jump of dispersion on the very early (prekinetic) stage could give a substantial contribution to 
overall dispersion. In this case the distribution function is governed by a Fokker-Planck equa- 
tion with time dependent diffusion coefficient. Solutions of the Fokker-Planck equation with 
quasilinear and time dependent diffusion coefficients were found numerically, and an analytical 
approximation for them was proposed as well. 

The numerical experiment gives the power law behavior of the dispersion at simulation 
times, and solutions of the Fokker-Planck equation shows it for extended interval. However the 



exponent is not unique for all spectra of the same type, as this follows from scaling consideration 
[4], but depends on the particular spectrum. 



2 Numerical model 

We consider the motion of noninteracting particles in an external random electric field. The 
potential of the field is taken as a superposition of M waves [5] 

M 
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with fixed frequency u> and wave numbers from the interval (k — 2.5 Ak, k + 2.5Ak). The total 
intensity of the field tp\ is distributed between the partial waves according to the Gaussian 

2 2 2 5k (ki-k \ 2 



where k® is the central wave number, Ak is the width of the spectrum, 5k = ki + i — ki = 5Ak/M. 

Each realization of the field (1) is characterized by unique set of random phases {ai}. For 
a given realization of the potential (1) the equations of particle motion, 

x — v, 

v = ^E(x, t), E(x, t) = -^<p(x, t), (3) 

with the initial conditions x(0) = 0, v(0) = vq are integrated numerically to obtain the particle 
trajectory in velocity space v(t). Average particle velocity v(t) and velocity dispersion (Av 2 ) t 
are found by averaging over realization 
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In simulations length is normalized to 27r//c , time to 2ty/uj, dimensionless potential and 
spectrum width are 

e k\ Ak 
a = — ip— and a = — — . 
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Kubo number, Q, which is the ratio of the correlation time to characteristic period of particle 
oscillations for this model is 

Q = V^/d- (5) 
Overlap parameter 
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much exceeds the unit for most harmonics j except of few at the wings of Gaussian distribution 
(2), and particle motion can be treated as stochastic. Note, that the random phase ensemble we 
used here does not provide stochastization by itself, but gives the explicit way for calculation 



of Euler correlation function of fields. According to Eqs. (1), (2) the correlation function for 
the potential ((f 2 ) x t is of the form 

(Akx) 2 

{¥ 2 )xt = V?o ex P 7 cos(^t - k x). (6) 

Obtained in simulation v(t) and {Av 2 ) t will be compared in Section 4 with numerical and 
approximate analytical solutions of the Fokker-Planck equation. 

3 Equation for distribution function and approximate 
analytical solution 

Introduce here the particle distribution function f(v,t) as microscopic distribution function 
averaged over random phase ensemble and integrated over the spatial variable x. As far as the 
averaging over the ensemble of random phase does not imply the averaging over any small but 
finite time scale, the distribution function is defined at all time scales, as well for t < r cor , i.e. 
in prekinetic stage. Generalized Fokker-Planck equation for f(v,t) in external fields could be 
obtained from Ref. 6 in the form 

with a time dependent diffusion coefficient 

2 * 

D(v,t)= (£) J(E 2 ) VTiT dr. (8) 
o 

For correlation function (6) it takes the form 

D(v, t) = ^ (j-(p Ak) J dr exp -(i Akvr) 2 

o 

x I ^1 + 2(^|) 2 - ^(Akv t) 2 ^ cos(cj - k v)r + 2k vr sin(c; - fc u)r|. (9) 

It will be shown that for moderate Kubo number the agreement with simulation is recovered 
by more accurate treatment of distribution function evolution on early stage t < r cor . Here the 
use of a time dependent diffusion coefficient is required. 

The asymptotic value of D(v,t) at large times t ^> r cor gives the well known quasilinear 
diffusion coefficient D qt (v) 

D ql (v) = D(v,t^oo). (10) 

When the correlation function is given by Eq. (6) the quasilinear diffusion coefficient takes the 
form 
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In the following section it will be shown that in the cases of narrow and/or high intensive 
spectrum it is important to retain the dependence of diffusion coefficient on time. 

Eqs. (7), (9) or (7), (11) which determine the evolution of f(v,t) with time dependent 
or, respectively, quasilinear diffusion coefficient are solved numerically. Initial condition for 



f(v,t = 0) was taken as a narrow Gaussian distribution that approximate S(v — v ). Then 
average velocity and dispersion are calculated as 



v (t) = / dvvf(v,t), 



(At; 2 ), 



dv(v-v(t)) 2 f(v,t). 
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They are compared in the following section with v(t) and (Av 2 )t obtained in simulation. 

In Fig. 1 the time dependent diffusion coefficient D(v,t) for narrow spectrum of external 
field along with its profiles are shown. 
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Fig. 1. Time dependent diffusion coefficient D(v,t) for times less than the correlation time 
T cor (left), sections of D(v,t) at t — 1, 3 and 100 (right), a = 0.01, d = 0.04. 



Time dependent diffusion coefficient (9) evolve from very broad distribution through oscillating 
regime to its asymptotic quasilinear value (11). 

Approximate WKB-type solution of the Fokker-Planck equation (7) with velocity dependent 
diffusion coefficient D q i{v) for the distribution function f(v,Vo,t) with initial conditions 
f(v, t = 0) = 5(v — vq) may be given in the form 

f(v,v ,t) = C(t)exp (-^pT), (14) 
with 

r v d,n 

Y(v,v )= / -^L= (15) 
and C(t) to be taken from the condition of normalization 
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This approximation was proposed in Ref . 4, however with other C (t) , which does not gives a 
proper time scaling of dispersion. 

For a time dependent diffusion coefficient the approximate WKB-type solution could be 
generalized as 

f(v, v , t) = C(t) exp (-\z\v, v , t)) , (17) 
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n (18) 
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Here, similarly to the previous case, C{t) should be defined from the condition of normalization. 

4 Comparison of simulation with numerical and analyt- 
ical solutions of the Fokker-Planck equation 

In this section results of simulations are compared with numerical and analytical solutions of 
the Fokker-Planck equation. For small Kubo numbers the solutions of Fokker-Planck equation 
with asymptotic quasilinear diffusion coefficient give the similar evolution of velocity dispersion 
and average velocity as the solutions with time dependent diffusion coefficient. In addition, for 
very small Kubo numbers at the beginning of evolution the velocity dispersion grows almost 
linearly. Whether Kubo number increase the deviation from the linear law due to dependence of 
diffusion coefficient on velocity becomes evident, however the solutions with D q i(v) and D(v,t) 
are still rather close. For a moderate Kubo numbers of the order of the unit, the difference 
between solutions with time dependent and asymptotic diffusion coefficient becomes noticeable. 

In Fig. 2 the curves obtained in simulation for a wide spectrum and moderate field are com- 
pared with the numerical solution of Fokker-Planck equation and the WKB solution (17), (18). 
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Fig. 2. Dispersion (Av 2 ) and average velocity v for a wide spectrum and moderate field, 
d = OA, a = 0.01 and Vo — 1. Kubo number Q = 0.25. Simulation and solution of 
Fokker-Planck equation with diffusion coefficient D(v,t) are compared with WKB solution 

(17), (18). 



In Fig. 3 is shown how WKB solution reproduces the early evolution of (Av 2 ) in the case of a 
narrow spectrum and moderate field with initial jump of dispersion. 

In the cases of broad spectrum and low intensity particles for long time slowly diffuse on 
small distance on v, which is less than halfwidth of D q i(v); and this time is enough for D(v,t) 
to evolve to its asymptotic value. In the opposite case, corresponding to Fig. 3, particles on 
times substantially less then r cor , while D(v,t) is a broad in v (c.f. Fig. 1), diffuse at large 
distance which is more than halfwidth of D q i(v). 
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Fig. 3. Dispersion (Av 2 ) for a narrow spectrum and moderate field, d = 0.04, a = 0.01, 
Vo = 1.0. Kubo number Q = 2.5. Simulation and solution of Fokker-Planck equation with 
D{y,t) are compared with WKB solution (17), (18) in a small time scale. 



The above examples were given for particles which initial velocity are not so far from the 
phase velocity of the central harmonic. Such particles in each instant are in resonance with 
some harmonic of considerable intensity and diffusion prevail over oscillations. For particles 
which initial velocities are far from resonance with intensive harmonics the diffusivity is small 
and oscillations become more distinct. In Fig. 4 dispersion are given for nonresonant particle. 
The curves obtained in simulation, as numerical solution of the Fokker-Planck equation with 
time dependent diffusion coefficient and in WKB approximation are shown. 
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Fig. 4- Dispersion of nonresonant particles (Av 2 ) for a narrow spectrum and weak field, 

d = 0.04, a = 0.0001, Vq = 1.2. Kubo number Q = 0.25. Simulation and solution of 
Fokker-Planck equation with diffusion coefficient D(v,t) is compared with WKB solution 

(17), (18) on a small time scales. 



5 Power law dispersion 

Simulation shows the velocity dispersion obeys a power law dependence on time. Numerical 
solutions of Fokker-Planck equation give the same power law, as simulation, and it easily could 
be calculated for much longer time. Such power law is also recovered from WKB approximation, 
and in this case it could be related to power law dependence of normalizing constants C(t). 
The example with Q = 0.079 is given in Ref. 7. Here, in Fig. 5 the plot is given for the different 



Kubo number, Q = 0.25, for numerical and WKB solutions of Fokker-Planck equation along 
with time dependence of G — C(t)t q (note that q = p/2). 
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Fig. 5. U = (Av 2 )/t p against t. Power law dispersion for a wide spectrum, p = 0.29 (left), and 
G — C(t) t q against t, q = p/2 = 0.145 (right); d = OA, a = 0.01, vq — 1, Kubo number 

Q = 0.25. 



The velocity dispersion shows power law time dependence, with exponent dependent on partic- 
ular spectrum. In terms of WKB solution it is related to power law dependence of normalizing 
constant C(t). 



Conclusions 

For small field intensity and wide spectrum (Kubo number less than the unit) the solution of the 
Fokker-Planck equation with quasilinear diffusion coefficient gives good agreement with results 
of numerical experiment. To have a consistency for high intensity and/or narrow spectrum 
(Kubo number of the order or larger than the unit) the generalization of the Fokker-Planck 
equation is to be done by introducing a time dependent diffusion coefficient. An analytical 
approximation for such solutions is proposed. Velocity dispersion manifests power law time 
dependence and the exponent is dependent on the spectrum. 
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